Scaling Between Periodic Anderson and Kondo Lattice Models 
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Continuous-Time Quantum Monte Carlo (CT-QMC) method combined with Dynamical Mean 
Field Theory (DMFT) is used to calculate both Periodic Anderson Model (PAM) and Kondo Lat- 
tice Model (KLM). Different parameter sets of both models are connected by the Schrieffer- Wolff 
transformation. For degeneracy N = 2, a, special particle-hole symmetric case of PAM at half filling 
which always fixes one electron per impurity site is compared with the results of the KLM. We 
find a good mapping between PAM and KLM in the limit of large on-site Hubbard interaction U 
for different properties like self-energy, quasiparticle residue and susceptibility. This allows us to 
extract quasiparticle mass renormalizations for the f electrons directly from KLM. The method is 
further applied to higher degenerate case and to realsitic heavy fermion system CeRhlns in which 
the estimate of the Sommerfeld coefficient is proven to be close to the experimental value. 
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I. INTRODUCTION 

Computational study of heavy fermion materials^^ is 
a challenging theoretical problem. These systems are 
a subset of intermetallic compounds that have a low- 
temperature specific heat whose hnear term is up to 
1000 times larger than the value expected from the free- 
electron theory. The heavy fermion behavior has been 
found in rare earth and actinide metal compounds at very 
low temperatures (typically less than 10 K) in a broad 
variety of states including metallic, superconducting, in- 
sulating and magnetic states^. 

The physics of the heavy fermion systems is controlled 
by the antiferromagnetic interactions of local magnetic 
moments residing on the rare earth or actinide atoms 
with the sea of conduction electrons. The theoretical 
problem of a localized spin interacting with the conduc- 
tion electrons is the celebrated Kondo problem'^"* whose 
solution is one of the outstanding achievements of many- 
body physics. It describes how the local spin is com- 
pensated as the temperature falls below a characteris- 
tic Kondo temperature. Something similar occurs in the 
heavy fermion materials which represents an array of such 
spins forming a Kondo lattice. 

In this regime, each f orbital is occupied by a fixed 
number of electrons, and all types of charge fluctuations 
are approximately frozen due to a large Coulomb repul- 
sion penalty that the system pays when the electron is 
added/removed from the shell. Therefore the low-energy 
degrees of freedom are provided by localized spins only 
and the corresponding model is known as the Kondo Lat- 
tice Model (KLMf-i'^. The KLM effective Hamihonian 
is obtained by using a second-order perturbation with 
respect to hybridizatiorti^ of a more general Periodic An- 
derson Model (PAM}i2ri^where the localized f-electrons 
can exchange with the conduction electrons bath thus al- 
lowing both charge and spin fluctuations to occur. The 



introduction of the limit of infinite dimensions and sub- 
sequent development of the dynamical mean field the- 
ory (DMFT) has allowed to study the properties of both 
models in a systematic manner—"— 

Due to developments in the electronic structure the- 
ory for strongly correlated systems based on a combi- 
nation of Density Functional Theory (DFT) in its local 
density approximation (LDA) and DMFT^-, studies of 
real heavy fermion materials have recently appeared in 
the literature^"—. Here the development of Continu- 
ous Time Quantum Monte Carlo Method (CT-QMC) for 
solving corresponding Anderson Impurity problem has 
played a central role^"— . These calculations are ex- 
tremely computationally demanding especially for the f 
elements such as Plutoniu m^^i^^ where a large number of 
atomic states needs to be kept in the calculation. 

We have recently proposed a simplified approaches 
where instead of full solution of the Anderson impu- 
rity model, a corresponding Kondo Impurity (or more 
general Coqblin-Schrieffer Impurity^) is studied to ex- 
plore low-energy physics of heavy fermion materials^i^ 
using most recently developed CT-QMC algorithm for 
this problem^"—. In this regard, an interesting ques- 
tion arises on how exactly the scaling between Ander- 
son and Kondo impurity models occurs and whether 
the low-energy properties of heavy fermion systems such 
as electronic mass enhancement and associated with it 
linear specific heat coefficient can be recovered from a 
restricted solution provided by the conduction electron 
self-energies available within the KLM. Such scaling be- 
havior has been explored^ for the temperature depen- 
dent susceptibility of the symmetric Anderson model us- 
ing numerical renormalization group techniques'", where 
a precise mapping has been found to spin-i Kondo 
Hamiltonian. Here we explore a similar mapping between 
singic-particle functions such as the self-energy where 
upon increasing the value of on-site Coulomb repulsion 
U, we report a convergence of the conduction electron 
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self-energy extracted from the solution of the PAM to the 
one obtained within KLM. We use CT-QMC method for 
the corresponding impurity models and the Dynamical 
Mean Field Theory for achieving self-consistent solution 
of the lattice problem. We utilize an inverse relationship 
to extract the f-electron self-energies and monitor how 
the low-frequency behavior of the AIM converges to the 
KLM limit. Our obtained mapping allows the extrac- 
tion of the mass renormalization of heavy quasiparticles 
directly from the solution of the Kondo lattice Hamilto- 
nian. 

As an illustration, we consider an electronic structure 
of CeRhlns where we compute hybridization functions of 
the f-electrons with conduction bath and evaluate Kondo 
exchange coupling. We subsequently solve the Kondo 
Lattice Model with CT-QMC and DMFT, compute con- 
duction electron self-energies, and then use the inverse 
mapping obtained from our analysis of the model Hamil- 
tonian to evaluate electronic mass enhancement and spe- 
cific heat coefficient of this system. Our theoretical re- 
sults are compared with available experimental data. 

This article is organized as follows. In Section II, we 
discuss the mapping between periodic Anderson model 
and Kondo lattice model in the limit of large U, and pro- 
vide the results for electronic self-energies, quasiparticle 
residues and susceptibilities. In Section III, application is 
presented to evaluate electronic mass enhancement and 
Sommerfeld's coefficient for CeRhlus. Section IV is the 
conclusion. 



Hamiltonian which eliminates the hybridization term 
in first order by Schrieffer- Wolff transformation. The 
second-order in V Hamiltonian is a famous Kondo Lat- 
tice Hamiltonian^"^'' that describes interaction between 
spins of localized and conduction electrons 

^^KLM = ^ <^\Lac\^„CVia + Jx'^Si ■ a (2) 
k(T i 

where Si represents the localized spin of the f electron at 
the i site, the a is the spin operator of the itinerant con- 
duction electron and Jk is the Kondo coupling constant. 



Jk = V^- 
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1 



-Ef ef + U' 



(3) 



For the half-filled case, e/ = — ?7/2, and Jk is simplified 
to W'^/U. 



B. Dynamical Mean Field Theory 

Solutions of both models in a general case represent a 
complicated numerical problem. Using dynamical mean 
field theory, the algorithm breaks down into (i) the solu- 
tion of the corresponding (Anderson or Kondo) impurity 
problem and (ii) the self-consistency loop over hybridiza- 
tion functions which enforces lattice periodicity'^^. 

For the periodic Anderson model, the DMFT evaluates 
the local Green function for heavy electrons 



II. MODEL CALCULATION 
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A. Periodic Anderson and Kondo Lattice Models 

One of the popular models to describe the physics 
of heavy fermion materials is a periodic Anderson 
mode&rM,, The effective Hamiltonian is given by 



k(T 



zk(T 



(1) 



where cj^^ (cko-) creates (destructs) a conduction electron 
with momentum k, spin (and orbital) a and dispersion 
fko-; fia ificr) crcatcs (destructs) an / electron with spin 
tj and energy e/ on site i; nj^ ('^{^) is the number operator 
for / electron at lattice site i with spin up (down); U is 
the on-site Coulomb repulsion; Vk is the hybridization 
between / electrons and conduction electrons which we 
assume to be k-independcnt, Vk — V, for simplicity. 

In systems where the Hubbard U is large, the charge 
fluctuations become effectively frozen and the ground 
state wave function has a little weight of configurations 
with the number of f-electrons different from its aver- 
age number fif. This results in transforming the PAM 



Then, the bath Green function is defined 

and used as an input to the impurity solver. The latter 
produces an impurity Green function 
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from where, a new self-energy can be found 

S/(za;„) = Go ^(iw„) - G^™'^''"^ (ia;„). 

The process is repeated by recalculating the lattice Green 
function with the new self-energy. The self-consistency 
condition is when 



G(""^)(zc.„)-Gp(zc.„). 

The Kondo lattice Hamiltonian can be obtained by 
considering the limit V"^— >oo,[/— >cx3,ey— >— cx) while 
keeping V'^/ef — const. First, define a local Green func- 
tion for conduction electrons 



G('°=)(i^„) = Y - £k - S,(ZC^„)]-1 
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where conduction electron self energy is given by 



N=2, T=0.01. J|,=0.3 



Sc(iw„) 



Thus, one can iterate over conduction electron quantities: 
The bath Green function 



serves as the input to the Kondo impurity solver. The 
latter produces an impurity Green function Gc*™^^ («cj„) 
from which the new conduction electron self-energy 
is found. The process is repeated by reevaluating 
Gc°'^''(ia;„).The self-consistency condition is when 

G(™f)(zw„) =G('°=)(zc<j„)- 

Several powerful methods such as exact diagonaliza- 
tion or numerical renormalization group techniques have 
been developed in the past to deal with the impurity 
models. In this work we utilize a Continuous Time Quan- 
tum Monte Carlo method^^"— that was originally pro- 
posed to deal with Anderson impurities but has been 
recently generalized for Kondo (Coqblin-Schrieffer) type 
of impurities^. 

The density of states of conduction electrons is an in- 
put to the simulation. Despite realistic materials may 
have complex band structures, we use a simple constant 
density of states to gain the physical insight from these 
calculations. The half-bandwidth D is set to 1 which 
provides the corresponding units. As we are looking for 
a mapping between the two models in the regime of large 
U, we first fix the Kondo coupling Jk to some predeter- 
mined value. There are typically two phases that emerge 
in the KLM: the antiferromagnetic RKKY phase and the 
paramagnetic Fermi liquid phase, which compete with 
each other on the scale of . We are mainly inter- 

ested in the Fermi liquid behavior and consider the value 
of Jk = 0.3 in all our calculations. Second, we study two 
cases with the effective f-electron degeneracies N = 2 
and iV = 4. For the case N = 2, the only non-trivial oc- 
cupancy of the f-orbital is 1 which for the particle-hole 
symmetric placement of the conduction electron band 
results in the condition ej = —U/2, for the / orbital 
to be half-filled. Although the system becomes a band 
(Kondo) insulator in this case, the f-electrons states are 
strongly renormalzed by correlations which is the basis 
for the comparison of these two models. 



C. DMFT Solutions for Af = 2 

1. Self- energies 

We first obtain the f-electron self-energy from the 
PAM calculation. Then we extract the conduction elec- 
tron self-energy and compare with the data of the KLM 
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FIG. 1: (color online) Conduction electron self-energy of the 
periodic Anderson model with N = 2 calculated for several 
values of U and the conduction electron self-energy of the 
Kondo lattice model that correspoinds to U oo limit. 



simulation. We present such comparison in Fig. [T] where 
we plot 9I]c(*<^n) for several values of U and U ^ oo 
limit corresponding to KLM. We monitor a slow conver- 
gence of the PAM self-energy towards its KLM value, 
although even for [/ = 10 the discrepancy between the 
two is still noticeable. 

We now address the question of ?7 — )■ oo limit for the 
f-electron self-energies, which numerically corresponds 
to the Kondo regime, and compare these data with our 
scaling behavior established analytically. First, notice 
that the low-frequency expansion for both f-electron and 
conduction electron quantities can be derived without a 
problem 

T.f,c{iuJn) = 5]/,c(0) -f iw„(l - Z^J) 
which leads us to 



S/(0) 
2/ 



F2 



Se(0) 



[5KS,(0)]2 



1/2 



1 - Zr, 



(4) 



This formula has been used in our recent LDA-f DMFT 
work to extract mass renormalization parameters in sev- 
eral heavy fermion compounds—- using simulations with 
the Kondo lattice. 

Unfortunately, this approach will not work for the 
model considered here since for the particle-hole sym- 
metric case of the Kondo insulator, the conduction elec- 
tron self-energies diverge to produce an energy gap in 
the excitational spectrum. We therefore look for a scal- 
ing behavior in a different way. We write 



(5) 
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N=2, T=0.01. J|,=0.3 
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FIG. 2: (color online) The convergence for the f-electron self- 
energy obtained from the periodic Anderson model with = 
2 upon increase in the interaction U. The limiting behavior of 
this quantity extracted from the solution of the Kondo lattice 
model is also plotted. 



and noticing that we target the [/ — > cxd limit, we replace 
V'^ with jJkU, and divide both parts by U. For the 
imaginary part we obtain the following scaling behavior 



U 



4|Sc(ia;„ 



(6) 



that expresses the large U limit of the PAM self-energy 
via the quantities available within KLM. 

The self-energies from both models can now be directly 
compared. Fig. [2]shows the behavior of SE/(iw„)/J7 for 
several values of U together with the corresponding data 
extracted from KLM. From the figure we see that as U 
increases the PAM self-energy converges to that of the 
KLM. The plot actually includes both the intermediate 
regime and the Kondo regime. When U < 2 the two 
Hubbard bands are within the conduction electron band 
which has a bandwidth of 2 in our units. In this case the 
self-energy deviates more than the data for t/ ^ 2. As 
U goes to 8, the self-energies from both models collapse 
at low energies. 

Finally, we notice here that the [/ — >■ oo limit pro- 
duces QY,f(iuJn)/U that grows linearly with the fre- 
quency which is exactly the case of the self-energies ob- 
tained in slave-boson type of methods for solving the 
impurity problematic or within quasiparticle Gutzwiller 
approximation^. 



2. Quasiparticle Residues 

Our self-energy results show qualitative convergence 
between the two models. To get quantitative agreement 
we further check the quasiparticle residues that we ex- 
tract from the low frequency behavior of the self-energies. 
These are related to renormalized effective masses for the 



FIG. 3: (color online) Comparison between quasiparticle 
residues Zf calculated as a function of Hubbard U using pe- 
riodic Anderson model with N — 2, Kondo lattice model as 
well as a slave-boson method described in Ref. 5a. 



quasiparticles responsible for the enhanced specific heat 
coefficient which is one of the central properties of sys- 
tems with heavy fermions. 

The CT-QMC algorithm works on an imaginary time 
axis, which after the Fourier transformation gives us the 
data on the imaginary frequency axis. Despite problems 
associated with numerical noise that prevents us to ex- 
tract accurate data at real frequencies using analytical 
continuation, we can find the quasiparticle residue from 
the imaginary axis data as follows: 



— T— |w„^oj • 



(7) 



For the KLM, according to Eq. ([5]), the Zf can be written 
as: 



KLM 
'1 



AttT |Ec(«7rr)|2 
' JkU 3I]c(j7rr) 



(8) 



Although this expression is actually valid only for U = 
cx), where zj^^ becomes zero, we expect that it gives 
an approximate value for PAM with U < oo. 

We present our comparisons between the two quanti- 
ties in Fig. [3] where we plot the quasiparticle residue 
extracted from PAM as a function of U as well as the 
one extracted from the KLM according to Eq. ([5]). Also 
plotted for comparison is the quasiparticle residue calcu- 
lated using the slave-boson method as described in Ref. 

m 

From the plots, we can quantitatively see the conver- 
gence from the result obtained from PAM to the one 
obtained by KLM. When U reaches 10, the KLM over- 
estimates Zf of the PAM data by about 30%. While 
the slave-boson method, a very fast calculation, demon- 
strates a similar behavior, it overestimates the PAM data 
by about 100%. 
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FIG. 4: (color online) Calculated using the periodic Anderson 
model with N = 2 temperature dependence of spin suscepti- 
bility (times the temperature) for several values of Hubbard 
U as well as the data extracted for the Kondo lattice model. 



FIG. 5: (color online) Conduction electron self-energy of the 
periodic Anderson model with A'^ = 4 calculated for several 
values of U and the conduction electron self-energy of the 
Kondo lattice model that correspoinds to U oo limit. 



3. Susceptibilities 

Since the KLM freezes the spatial fluctuations but 
keeps the essential magnetic properties, the correspond- 
ing spin susceptibilities x(T) should agree between the 
PAM and KLM at least for low temperatures. In- 
deed, these quantities have been calculated and compared 
against each other for the corresponding impurity mod- 
els using numerical renormalization group methods long 
time agc**^ where a precise mapping between Anderson 
impurity and spin-i Kondo impurity has been observed. 

We present our own calculations in Fig. 2] where we 
plot Tx{T) against the temperature for several values of 
U calculated using the PAM as well as the data extracted 
from the KLM. The comparison shows a nice conver- 
gence for susceptibility within our chosen temperature 
range. The deviation may result from the combination 
of thermal effect and charge fluctuations. We see that the 
convergence is worse here than that for the quasiparticle 
residues discussed earlier, but we believe that the sus- 
ceptibilities should map precisely to each other at lower 
temperatures. 



D. DMFT Solutions for A = 4 

The N = 2 case with one localized electron leads to 
the Kondo insulator state which is topologically special. 
However, our method can be generalized to larger orbital 
degeneracy, where any integer occupancy of the f-shell 
can be explored. Also the Coqblin-Schrieffer model is 
more favorable. Below we consider the case with = 4 
and Hf — I which is away from particle-hole symme- 
try. As is the case with A^ = 2, we fix the value of the 
Kondo coupling Jk to 0.3. For each value of U that 



we input to the PAM calculation, there are two remain- 
ing parameters, the impurity level e/, and the value of 
hybridization that should be searched for to obtain 
Uf — 1, Jk = 0.3. 

Conduction electron self-energies 3S]c(»w„) calculated 
within PAM for several values of U as well as within KLM 
corresponding to U — > oo limit are compared in FiglS] 
We see that the convergency of $5I]c(iw„) is rather slow 
when U increases similar to the N = 2 case presented in 
Fig. [TJ The low frequency behavior of 3Sc(ia;„) shows 
that the hybridization gap is no longer opened at the 
Fermi energy and the system remains metallic contrary 
to the particle-hole symmetric case of the Kondo insula- 
tor where 9Sc(?Wn) diverges at za;„ — > as seen on Fig. 
[T] Here, the low-frequency slopes determine quasipar- 
ticle residues Zc for conduction electrons which display 
a somewhat faster convergence to the Kondo limit upon 
increase in U. 

To compare how Sy(ia;„) scales to the Kondo limit, we 
start from Eq. ([5]) , take its imaginary part and divide 
by U on both sides. Using Eq. ([3]), the formula becomes: 



u 



+ JKTlil - ri) 



|Sc(iw„) 



(9) 



- is a dimensionless parameter which affects 



where rj — jj 
the electron counting at impurity site. 

Fig. (jni presents the behavior of 3I]/(iaj„)/C/ calcu- 
lated within PAM for several values of U together with 
the corresponding data extracted from KLM. From the 
figure we see that as U increases, the PAM self-energy 
maps into its U — > oo limit of the KLM. 

We have finally extracted the values for quasiparti- 
cle residues from the low-frequency slopes of 3Ej:(iw„) 
which can be compared with the values of that we ob- 
tain from the KLM calculation either using the approach 
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N=4, T=0.01, J^=0.3 
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FIG. 6: (color online) The convergence for the f-electron self- 
energy obtained from the periodic Anderson model with A'^ = 
4 upon increase in the interaction U. The limiting behavior of 
this quantity extracted from the solution of the Kondo lattice 
model is also plotted. 
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FIG. 7: (color online) Comparison between quasiparticle 
residues Zf calculated as a function of Hubbard U using pe- 
riodic Anderson model with A" = 4, and the values extracted 
from the Kondo lattice model using two different approaches 
described in text. 



that leads us to Eq.([8]) or using the low-frequency be- 
havior of '^T,c{iujn) that leads us to Eq. (j4]). We present 
such comparison in Fig. [7] where the behavior of Zf is 
plotted against Hubbard U. We see that starting from 
U = 10, the quasiparticle residues computed from PAAI 
and KLM become very close to each other. 



III. APPLICATION TO CERHIN5 

Realistic heavy fermion materials have much more 
complicated electronic structures than we used in our 
model calculations. The f-orbitals are 14-fold degen- 



TABLE I: Calculated LDA density of states at the Fermi en- 
ergy, N{0) (states/eV/cell), estamted quasiparticle residue 
Zf for the f-electrons as well as predicted and experimen- 
tal values of the Sommerfeld coefBcient 7 (mJ/mol/K^) for 
CeRhlns. 



materials 7V(0)lda zj 7 'yexp 
CeRhlug 2.21 0.01255 414 400° 

^RifU 

erate and split in the presence of spin-orbit coupling 
and crystal-field effects. This makes LDA-hDMFT 
calculations with full solution of Anderson impurity 
problem dramatically heavy. With the advantage of 
reaching lower temperature range, LDA-I-DMFT sim- 
ulations with Kondo (or Coqblin-SchriefFcr) impurity 
was established^^. In this method, the first step is to 
find the a local hybridization function A.a{e) between 
fully localized / and conduction spd electrons, where 
subscript a refers to a particular representation that 
tries to diagonalize a general matrix AmVm(T(e) by tak- 
ing advantage of spin-orbit and crystal-field symme- 
tries. This is achieved by using the so called Hub- 
bard I approximatioEk^ where purely atomic f-electron 
self-energy is inputted to the LDA-I-DMFT calculation. 
Knowing Aa(e), the Kondo coupling constant Jk and 
the initial Green's function of conduction electrons can 
be extracted. 



Ga (iuJn) 



-Dcutotf C>A /"-Dcutotf 

de^^==^/ / de3A„(e), 



and 



T/2 ^ 



1 f""^"'"" , 5Aa(e) 
ae- 



(10) 

(11) 



These are required for solving realistic Kondo lattice 
problem with Jk calculated from the Schrieffer- Wolff 
transformation, Eq. ([3|). 

Our target material is CeRhlns which is believed to 
have most localized f-electrons in the 115 family. The 
spin-orbit coupling and crystal fields of the tetragonal 
structure effectively reduce the degeneracy and make the 
r7 doublet of the j = 5/2 state to be the ground state. 
Therefore at very low temperatures, a single localized f- 
electron resides at the Fy doublet and we have the N = 2 
case discussed above in the model calculation. 

A general LDA-HDMFT calculation for this material 
has been done in the former work'^^ with ef — — 2.5eV 
and U = 5eV which are the typical values for Ce-based 
compounds. Here we provide an analysis of its low en- 
ergy physical properties. Fig. [5]shows calculated conduc- 
tion electron self-energies, $HEc(e) and 5Sc(e), for the T^ 
state from our LDA-I-DMFT simulation with the Kondo 
lattice. We clearly see that the imaginary part of Sc 
tends to diverge at least to T p» 23K although this be- 
havior could result from the temperature being not low 
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FIG. 8: (color online) Conduction electron self-energies 
$HEc(e) (top plot) and SfSc(e) (bottom plot) of Ft states for 
CeRhlns at temperature T = 0.002 eV(fa23K) calculated 
using the LDA-f DMFT method with the Kondo lattice. 



enough in our simulation. However, a conventional low- 



frequency expansion of the self-energy and the connec- 
tion between Zf and Zc, Eq. ([4]), cannot be utilized to 
estimate the quasiparticle residue for the f electrons and 
the Sommerfeld coefficient 7. 

Here we use our mapping method to extract Zf exactly 
as we illustrated for our model calculation. In this way, 
we first estimate the quasiparticle residue Zf and, second, 
evaluate the renormalized density of states at the Fermi 
level A^(0)ofr = A^(0)lda/^/- Then the Sommerfeld co- 
efficient can be found 



7 = -7riVeff(0). 



(12) 



All calculated properties are summarized in Table HI It 
can be seen that our estimate for 7 is very close to the 
experimental value which indicates that our simulation 
is sufficiently accurate to describe this material. 

IV. CONCLUSION 

We have studied a mapping of periodic Anderson 
model to the Kondo lattice model in the limit of t/ — !■ 00 
for single-particle functions such as the self-energy. The 
crossover occurs at the values of interaction U = \QD 
where the models become equivalent. This allowed us to 
map the quasiparticle residue zj ol the f-electrons and 
extract its values directly from the Kondo lattice model. 
We applied the method to realistic heavy fermion system 
CeRhlns where our estimates for the Sommerfeld coeffi- 
cient agree well with the experiment. 
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